General relations for quantum gases in two and three dimensions. 

II. Bosons and mixtures 
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We derive exact general relations between various observables for N bosons with zero-range in- 
teractions, in two or three dimensions, in an arbitrary external potential. Some of our results are 
analogous to relations derived previously for two-component fermions, and involve derivatives of the 
energy with respect to the two-body s-wave scattering length a. Moreover, in the three-dimensional 
case, where the Efimov effect takes place, the interactions are characterized not only by a, but 
also by a three-body parameter Rt. We then find additional relations which involve the derivative 
of the energy with respect Rt. In short, this derivative gives the probability to find three parti- 
cles close to each other. Although it is evaluated for a totally loss-less model, it remarkably also 
gives the three-body loss rate always present in experiments (due to three-body recombination to 
deeply bound diatomic molecules), at least in the limit where the so-called inelasticity parameter 
r\ is small enough. As an application, we obtain, within the zero-range model and to first order in 
7/, an analytic expression for the three-body loss rate constant for a non-degenerate Bose gas with 
infinite scattering length. We also discuss the generalization to arbitrary mixtures of bosons and/or 
fermions. 

PACS numbers: 67.85.-d 



I. INTRODUCTION 

Ultracold gases with resonant interactions, that is hav- 
ing a s-wave scattering length much larger in absolute 
value than the interaction range, can now be studied 
experimentally thanks to the broad magnetic Feshbach 
resonances, not only with fermionic atoms [l[ but also 
with bosonic atoms or mixtures @. In this reso- 

nant regime, one can neglect the range of the interaction, 
which is equivalent to replacing the interaction with con- 
tact conditions on the A^-body wavefunction: In 3D, this 
constitutes the so-called zero-range model 

ill/ 

that 

can also be defined in 2D (see e.g. |13l - ilq |). and of course 
in ID [T3, EH- I n eacn dimension, these models include 
a length, the so-called d-dimensional scattering length a. 
In three dimensions, when the Efimov effect occurs @, 
an additional length has to be introduced, the so-called 
three-body parameter [l9j . 

For the zero-range models, it was gradually realized 
that several observables, such as the short distance be- 
havior of the pair distribution function g^(r) or the 
tail of the momentum distribution n(k), can be re- 
lated to derivatives of the energy with respect to the 
d-dimensional scattering length a. In ID, the value of 
g( 2 \0) was directly related to such a derivative by the 
Hellmann-Feynman theorem [l7j]; the coefficient of the 
leading 1/k 4 term in n(k) at large k was then related to 
the singular behavior of the wavefunction for two close 
particles, and ultimately to g^(0), by general properties 
of the Fourier transform [2(|. In 3D, for spin-1/2 fermions 
(where the Efimov effect does not occur), an extension of 
the ID relations was obtained by a variety of techniques 
[21I - I26I ] . including the original ID techniques. General- 



izations were then obtained for 2D systems, for fermions 
or bosons [IMoj]. 

This is the second of a series of two articles on such 
general relations. The first one covered two-component 
fermions (Ref. (30[, hereafter referred to as Article I). 
Here, we consider single-component bosons, as well as 
mixtures. In the 3D case, remarkably, the Efimov ef- 
fect leads to modifications or even breakdown of some 
relations, and to the appearance of additional relations 
involving the derivative of the energy with respect to 
the three-body parameter R t . Several of the results pre- 
sented here were already contained in [3l| and rederived 
in [32| with a different technique, that allowed the au- 
thors of [32f to obtain still other Efimovian relations for 
N bosons [6a ]. 

The article is organized as follows. Section [XT] intro- 
duces the zero-range model and associated notations for 
the single-component bosons. Section IIIII presents rela- 
tions which are analogous to the fermionic ones. Addi- 
tional relations resulting from the Efimov effect are de- 
rived in Section IIVI As an application, the three-body 
loss rate of a non-degenerate Bose gas for an infinite scat- 
tering length is calculated in Section |Vj Finally the case 
of an arbitrary mixture is addressed in Section IVII For 
convenience, the obtained relations are displayed in Ta- 
bles I, II, III. 



II. MODEL AND NOTATIONS 

In 3D, the zero-range model imposes the Wigner- 
Bethe-Peierls contact condition on the A^-body wavefunc- 
tion: For any pair of particles when one takes the 



2 



limit of a vanishing distance = jr, — rj\ with a fixed the differential solid angles in 3D, we obtain |33l. 134 
value of the center of mass c.^ = (r, + ij-)/2 different 
from the positions of the other TV — 2 particles, the 
wavefunction has to behave as 



V>(ri, ...,r N )= ( — - - J Aijicij, (r k )k^i,j) + 0{nj) 
\ r ij a J 

(1) 

where a is the 3D scattering length. The a priori un- 
known functions Ay are determined from the fact that 
■0 solves the free Schrodinger's equation over the domain 
where the positions of the particles are two by two dis- 
tinct: Eip = Hip with 



N 



h 2 

2m 



(2) 



and U is the external potential. Also ip is normalized to 
unity. 

If there are three bosons or more, the Efimov effect 
occurs @, and the zero-range model has to be supple- 
mented by a three-body contact condition that involves 
a positive length, the three-body parameter R t : In the 
limit where three particles approach each other (that one 
can take to be particles I, 2 and 3 due to the bosonic 
symmetry), there exists a function B, hereafter called 
three-body regular part, such that 



•,rjv) ~ 9(T U r 3 ,T 3 )B(C,u,...,r N ) (3) 

R-yO 



where C = (ri + r 2 + ra)/3 is the center of mass of parti- 
cles 1,2 and 3, $ is the zero-energy three-body scattering 
state 



$(ri,r 2 ,r 3 ) 



1 

R 2 



sin 



I so I In 



,(«), (4) 



and where R, f2 are the hyperradius and the hyperangles 
associated with particles 1, 2 and 3. We take the limit 
R I — > in © for fixed fl and C (in analogy with the 
two-body contact condition). 

We recall the definition of R and £1: From the Ja- 
cobi coordinates r = r 2 — r± and p = (2r 3 — ri — r 2 ) /V3, 
one forms the six-component vector R = (r, p)/y/2; then, 
the hyperradius R — \/ (r 2 + p 2 )/2 is the norm of R, and 
f2 = R/i? is its direction that can be parametrized by five 
hyperangles, so that d 6 R = R 5 dRd 5 Q. In Eq. (gj), s = 
%■ 1.00623782510 ... is Efimov's transcendental number, it 
is the imaginary solution (with positive imaginary part) 
of scos(s7r/2) = (8/V3) sin(s7r/6); (j) So (fl) is the hyper- 
angular part of the Efimov trimers wavefunctions @, 
which, in the present case (single-component bosons), is 
given by so (f2) =AA(1 + Q)sinh[|s | (f - a)] /sin(2a) 
where Q = P13 + P23 and Py exchanges particles i 
and j, and where a = arctan(r/p). Here we in- 
troduced, for later convenience, a normalization fac- 
tor such that f d 5 n\(j) So (n)\ 2 = 1. Using / d 5 ft = 

J^ 2 da sin 2 a cos 2 a J d 2 f Jd 2 p, where d 2 f and d 2 p are 



67T 2 r 
N~ 2 = -. — r sinh(|sok/2) cosh(|s |7r/2) 

7T 4l7T 

+ |s |-sinh(|s K/2)- —7 =cosh(|s k/6) . (5) 
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For N — 3 particles, it is well established that this 
model is self-adjoint and that it is the zero-range limit of 
finite-range models, see e.g. [l2j and references therein. 
The fact that the zero-range (i.e. low-energy) regime can 
be described using the scattering length and a three-body 
parameter only is known as universality [TT| . For N = 
4, an accurate numerical study [35[ has shown, as was 
suggested by earlier ones |36l - [38| , that there is no need to 
introduce a four-body parameter in the zero-range limit, 
implying that the here considered zero-range model is 
self-adjoint for N — 4. Physically, this is related to the 
fact that the introduction of Rt, imposed by the three- 
body Efimov effect, breaks the separability of the 4-body 
problem at infinite scattering length and prevents the 
occurrence of a four-body Efimov effect for bosons [39| . 
Here we consider an arbitrary value of N such that the 
model is self-adjoint. 

In 2D, the zero-range model is a direct generalization 
of the 3D one, since one simply replaces the 3D zero- 
energy two-body scattering wavefunction — aT 1 by 
the 2D one ln(r,j/a), where a is now the 2D scattering 
length. For any pair of particles i and j, in the limit 
rij = I r i — r j | — > with = (r< + r,-)/2 fixed, the 
iV-body wavefunction satisfies in 2D: 

tp (ri, ...,r N ) = ln(r y /a)Ay(cy, (r k ) k ^ i}j ) + 0(nj). 

(6) 

There is no Efimov effect in 2D so that no additional 
parameter is required (40l - |42| . The Hamiltonian is the 
corresponding 2D version of ([2]). 



III. RELATIONS WHICH ARE ANALOGOUS 
TO THE FERMIONIC CASE 

A first set of relations is given in Table HI These 
relations and derivations are largely analogous to the 
fermionic case (which was treated in Article I). An 
obvious difference with the fermionic case is that there 
are no more spin indices in the pair distribution function 
<?( 2 ) and in the momentum distribution n(k). Accord- 
ingly we now have g( 2 \r,r') = (ip^ (r)^ (r')tp(r')^p(r)) = 
J d d n . . . d d r N |V(ri, . . • , r,v)| 2 . 8 (r - r<) S (r' - r,-), 

where ip is the bosonic field operator, and the momentum 
distribution is normalized as J n(k)d d k / (2ir) d = N. 
Apart from numerical prefactors, there are two more 
important differences which appear in the 3D case due 
to the Efimov effect. 

The first important difference is that the derivatives 
with respect to 1/a in [Tab. HI Eqs. (4a, 6a)] have to be 
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TABLE I: For single-component bosons, relations which are analogous to the fermionic case. In three dimensions, the derivatives 
are taken for a fixed three-body parameter Rt. As discussed in the text, in three dimensions, the relation between energy and 
momentum distribution is valid if the large cut-off limit A — > +oo exists, which is not the case for Efimovian states (i.e. 
eigenstates whose energy depends on R t ). The notation (A, A) is defined in Eq. {5J. 7 = 0.577215 ... is Euler's constant. 



taken for a fixed three-body parameter R t . This comes 
from the relation 



dE 



(A A), 



(7) 



with the notation (given for generality in dimension d): 



(A, A) 



E 

i<j 



J ( J| d d r k ) J d d c l AA l0 {c l0 ,{r k ) k ^ l 



_ (8) 
Eq. ([7]) was already obtained in |l2[ in the case N = 3. 
A simple way to derive it for any N is to use a cubic 
lattice model, of lattice spacing b, with purely on-site 
interactions characterized by a coupling constant go [see 
the Hamiltonian in Eq. (|14[) with ho = 0], adjusted to 
reproduce the correct scattering length |43|: 



go Ai:h 2 a 



d 3 k 



D 



(2tt) 3 H 2 k 2 



(9) 



where the wavevector k of a single particle plane wave 
on the lattice is restricted as usual to the first Brillouin 
zone D = (— j, ^) 3 . One then follows the same reasoning 
as in (Article I, Section V, Subsections C-D-E). The key 
point here is that, in the limit of b -C \a\, the three-body 
parameter corresponding to the lattice model is equal 
to a numerical constant times b [67|. Thus, varying the 
coupling constant go while keeping b fixed is equivalent 
to varying a while keeping R t fixed, so that 



dE _ ( dE \ d(-l/a) 
d 9a \d(-l/a)) Rt dg a 



(10) 



The left-hand side of (jTU)) is given by the Hellmann- 



Feynman theorem: 



J Tp 1 

dg 2 ^ 



N(N — 1) 



J2 6 3(iV - 1) |V(r,r,r 3 ,...,r J v)| 2 (11) 



where ip is the eigenstate wavefunction on the lattice. In 
the zero-range limit b <C \a\, ip has to match the contact 
condition ([T]): Its two-body regular part A12, defined as 



V>(r, r, r 3 , . . . , r N ) = (p(0)A 12 (r, r 3 , . . . , rjv) 



(12) 



with the correctly normalized zero-energy two-body lat- 
tice scattering wavefunction <f>(r) [<f>(r) = r^ 1 — a -1 +o(l) 
at r 3> b] , has to converge to the zero-range model regular 
part. Similarly, in the right-hand side of (jTU)) . the lattice 
model's (dE/d{— l/a))n t tends to the zero-range model's 
one if one takes the zero-range limit while keeping R t 
fixed [6a |. Finally, the last factor of (|10p can be evalu- 
ated from ©. Using the relation (j)(0) = — 4ttH 2 / (mgo) 
established in [3(J, we obtain [Tab.HlEq. (4a)]. The same 
lattice model reasoning explains why the second-order 
derivative in [Tab. HI Eq. (6a)] also has to be taken for a 
fixed R t . 

The second important difference with respect to the 
fermionic case is that the relation [Tab. HI Eq. (5a)] breaks 
down in general, and only holds for special states for 
which the infinite-cutoff limit A — > 00 exists (such as the 
universal states for 3 trapped bosons of [3, HU). This 
was overlooked in [28| . and was shown for an Efimov 
trimer in [34| . The correct relation valid for any AT- body 
state in presence of the Efimov effect was obtained in [32| ■ 
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IV. ADDITIONAL RELATIONS COMING 
FROM THE EFIMOV EFFECT 

In addition to modifying relations which already ex- 
isted for fermions, the Efimov effect gives rise to addi- 
tional relations, involving the derivative of the energy 
with respect to the logarithm of the three-body parame- 
ter. These relations are displayed in Table [TTl 



wave annihilation operator obey as usual [ck , c^,] = 
(27r) 3 5(k — k') provided that k and k' are restricted to 
the first Brillouin zone D. 

We then define the zero-energy three-body scattering 
state o (ri,r 2 ,r 3 ) as the solution of 22i a tt|<A)) = for 
a = oo, with the boundary condition 

<Mri,r 2 ,r 3 ) ~ $(n,r 2 ,r 3 ) (15) 



A. Derivative of the energy with respect to the 
three-body parameter 

Our first additional relation [Tab, fill Eq. (1)] expresses 
the derivative of the energy with respect to the three- 
body parameter R t in terms of the three-body regular 
part defined in Eq. This is similar to the relation 
([7]) between the derivative with respect to the scattering 
length and the (two-body) regular part [69| . We will first 
derive this relation using the zero-range model in the case 
TV = 3, and then using a lattice model for any N. 

1. Derivation using the zero-range model for three particles 

We consider two wavefunctions ipi, ip2, satisfying the 
two-body boundary condition (|TJ) with the same scat- 
tering length a, and the three-body boundary condition 
PI4[) with different three-body parameters R t \ 1 Rti- The 
corresponding three-body regular parts arc denoted by 
B%, £>2- We show in the Appendix lAl that 



{iI>x,H^)-{H^i,tP 2 ) 



h 2 3V3| 



so 



i 1 1 

Mln — 

*H1 



J d 3 C5nC)B 2 (C), (13) 



which yields [Tab. II, Eq. (1)] by choosing ifji as an eigen- 
state of energy Ei and taking the limit Rti — ¥ Rti [7Cj . 



2. Derivation using a lattice model 



in the limit where all interparticle distances tend to infin- 
ity. Here $ is the zero-range model's zero-energy scatter- 
ing state, given in Eq. (Q}. This defines the three-body 
parameter R t (b,ho) for the lattice model (since $ de- 
pends on R t ). The Hellmann-Feynman theorem writes: 



dE 



— =£y<(V>W^)(r)> 



N(N-l)(N-2) & 3(JV -%(r,r,r,r 4) 



r.v . l- s 



(16) 



For the lattice model we define the three-body regular 
part B through: 

vj(y, r, r, r 4 , . . . , rjv) = O (O, 0, 0) B(r, r 4 , . . . , r N ); (17) 

in the zero-range limit, we expect that this lattice model's 
regular part tends to the regular part of the zero-range 
model defined in Eqs. (|3I4|) . Thus, in the zero-range limit: 



dE 



f)h 

= iV(iV-l)(JV-2)|0 o (O,O,O)| 2 ' 



d(lnR t )J a " " "™ ' ' " KdQnRtj;,, 
x J d 3 rd 3 r 4 ...d 3 r N \B{r,r A ,...,r N )\ 2 . (18) 



It remains to evaluate the derivative of ho with respect 
to R t : This is achieved by applying (fT5j) to the case of 
an Efimov trimer in free space, where the regular part 
can be deduced from the known expression [34| for the 
normalized wavefunction. This yields [Tab. II, Eq. (1)]. 



We now derive [Tab. II, Eq. (1)] for all N using as in 
Sec. IIIII a cubic lattice model, except that the Hamilto- 
nian now contains a three-body interaction term (of cou- 
pling constant ho) allowing one to adjust the three-body 
parameter R t without changing the lattice spacing: 

+ ^(V^WXr) +h Q J2 fo 3 (V' t V' t ^ t ^)(r). 

r r 

(14) 

Here the bosonic field operator obeys discrete commu- 
tation relations [i>(r), ifj'(r')] = S rr > /b 3 and the plane 



B. Short-distance triplet distribution function 

Similarly to the pair distribution function 
g^ 2 \ one defines the triplet distribution function 
5 ( 3 )(ri,r 2 ,r 3 ) = (^t^)^^)^^)^^)^^)^^)), 
which is given in first quantization by N(N — 1)(N — 
2) Jd 3 r 4 ■ ..d 3 r N |i/-(n, . . . ,rAr)| 2 . In the limit R -> 
where the three positions ri,r 2 ,r 3 approach each 
other, the many-body wavefunction behaves according 
to ©. The result [Tab. HI Eq. (2)] then directly 
follows, using [Tab. UH Eq. (1)]. As a consequence, in 
a measurement of the positions of all the particles, the 
mean number of triplets of particles having a small 
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/ dE 

\d(\nR t ) 


\ = ^l s ol 2 jv(7V _ 1)( ^_ 2) [ d 3 C d 3 r 4 ...d 3 r N \B(C,r 4 ,...,r N )\ 2 (1) 


J 


^^» (ri , r2 , r3) ^ i $(riir2)r3 )r ( a(1 f flt) ) a ^; so|2 1 ( 2) 


„->o \d(\nR t )J a \ ao \ W 



TABLE II: For single-component bosons in 3D, additional relations coming from the Efimov effect. B is the three-body regular 
part of the iV-body wavef unction; is the triplet distribution function; T is the decay rate due to three-body losses and r\ is 
the corresponding inelasticity parameter (see text). 



hyperradius R is given by 



AWets(# < e) = T\ I d 3 r 1 d 3 r 3 d 3 r 3 g^{r 1 ,r 2 ,r 3 ) 



>o 2h 2 \s Q \ 2 \d(lnR. 



dE 



1 - Re 



(e/Rt) 2lM 
l + i\s \ 



(19) 



where we used the Jacobian ^^^p^ = 3\/3 and the 
division by 3! takes into account the indistinguishability 
of the particles within a triplet. 



C. Decay rate due to three-body losses 

In experiments, the cold atomic gases are only 
metastable: There exist deeply bound dimer states, that 
is with a binding energy of order h 2 /{mb 2 ), where b is 
the van der Waals length of the real atomic interaction. 
These deeply bound states can be populated by three- 
body collisions, which are strongly exothermic (with re- 
spect to the trapping potential depth) and thus lead to a 
net loss of atoms. Usually, one expects that these deeply 
bound dimer states have a vanishing small effect on the 
metastable many-body states for b — > 0; the metastable 
states then converge to stationary states described by the 
zero-range model. 

In presence of the Efimov effect, however, the probabil- 
ity Pdosc to find three particles within a distance b (e.g., 
with an hyperradius R < b) vanishes only as b 2 accord- 
ing to Eqs. (I3I4I19P - As the three-body loss rate scales 
as p r .]n se.h/mb 2 , it does not vanish in the zero-range limit 
0, Eof - Fortunately, one can still in that limit simply in- 
clude the losses by modifying the three-body boundary 
conditions [13, 51]: One keeps Eq. ^ with a modified $ 
deduced from Eq. (|4]) by the substitution 



sm 



|s |ln 



1 

2i 



e -r, e i\s \\n{R/R t ) 



. e ri e -i\sa\HR/Rt) 



(20) 



The so-called inelasticity parameter r\ > determines to 
which extent the reflection of the incoming hyperradial 
wave exp[— %\sq\ bx(R/Rt)] on the point R = (where the 



model-dependent short range three-body physics takes 
place) is elastic. 

In this work, we have considered so far the ideal case 
where r\ is strictly zero. We now show that this allows 
to access the decay rate due to three-body losses to first 
order in r\ by taking simply a derivative of the loss-less 
eigenenergies E. In a first approach, we simply assume 
that E^lni?^ is an analytic function of \nR t . As the sub- 
stitution (|20[) simply amounts to performing the change 



In R f — > In R f 



17] 

R 



(21) 



we conclude that the resulting eigenenergy for non-zero 
r) acquires an imaginary part —ihT/2 given to first order 
in t] by [Tab. HH Eq. (3)]. Furthermore, we have devel- 
oped an alternative approach, that relates for arbitrary 
r] the decay rate F to the integral of |i?| 2 , where B is 
defined by Eq. Q , see Appendix |Bj Combining this with 
[Tab. HH Eq. (1)] in the limit r\ — > reproduces the rela- 
tion [Tab. El Eq. (3)]. 



V. APPLICATION: THREE-BODY LOSS RATE 
CONSTANT FOR A BOSE GAS 

We consider a 3D Bose gas, in a cubic quantization box 
of volume V , at thermal equilibrium in the grand canon- 
ical ensemble and in the thermodynamic limit. Within 
the zero-range model, with a truncation of the three-body 
energy spectrum (that is introducing a lower energy cut- 
off, as discussed below), relation [Tab. HH Eq. (3)] can 
be used to obtain, to first order in the inelasticity pa- 
rameter 77, the three-body loss-constant L 3 customarily 
defined by 



— N = —L 3 nN 
at 



(22) 



where N is the mean particle number and n = N/V the 
mean density. Applying [Tab. [TTl Eq. (3)] to each many- 
body eigenstate, taking a truncated thermal average and 
keeping in mind that each loss event eliminates three par- 
ticles out of the system |7l| . we obtain 



-j—W = °) = 
drj 



G 



h\s Q \n 2 N \dQnR, 



(23) 
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where the derivative of the grand potential fl is taken for 
fixed chemical potential \i and temperature T [72j . 

To obtain analytical results, we restrict to the non- 
degenerate limit fi — > — oo, where the density vanishes, 
n\ 3 -> 0, with A = [2-k% 2 / {mk B T)] x / 2 the thermal de 
Broglie wavelength. One then can use the virial expan- 
sion 149 531: 



n{jt,T) 



V 
A 3 "' 



(24) 



with /3 — l/(fcsT), and b q only depends on g-body 
physics and temperature. The leading order contribu- 
tion that involves lni? t is thus for q = 3, so that 



12tt h\ 4 



db* 



— ( =0) -> 

d?7 nA 3 ->o | s 1 m \d(\nR t ) / T 



(25) 



where we used nA 3 ~ exp(/3/i). 

The coefficient b q can be deduced from the solution of 
the g-body problem. We thus restrict to the resonant 
case 1/a = 0, where the solution for q — 3 is known in 
free space ly . Due to separability in hyperspherical co- 
ordinates [5J] the solution is also known for the isotropic 
harmonic trap case [iil . l45j, w hich allows us to use the 
technique developed in [53l |55| to calculate 63: 



63 = 3 3/2 lim 



Z 3 y , 1 7 2 



(26) 



where Z q (uS) is the canonical partition function at tem- 
perature T for the system of q interacting bosons in the 



harmonic trapping potential U (r) 



j 2 r 2 . Since the 



center-of-mass is separable, Z3/Z1 simply equals the par- 
tition function Z 3 nt of the internal variables. The internal 
3-body eigenspectrum in the trap involves fully univer- 
sal states (not depending on R t ), and a single Efimovian 
channel with i?t-dependent eigenenergies E n (uj), n e Z, 
solving a transcendental equation. Within the boundary 
conditions (|3I4|1 . the sequence E n {ui) is unbounded be- 
low. To give a mathematical existence to thermal equilib- 
rium, we thus truncate the sequence, labelling the ground 
three-body state with the quantum number n = and 
then keeping only n > in the thermal average f73j . 
In the free space limit lj — > 0, this corresponds to a 
purely geometric spectrum of trimer states with a ratio 
exp(— 27r / 1 so I ) and a ground state Efimov trimer energy: 



Eo(w) 



2h 2 
>o mR 2 t 



(27) 



Given Et, this uniquely determines the three-body pa- 
rameter R t [z3|. This finally leads to 



db a 



dQnRt)J T 



33/2 

lim > ( 

k B T ^0 4-, 

n>0 



-0E n (u) dE n (uj) 
9(ln R t )' 



(28) 

Details of the calculation of that limit are exposed in 
Appendix [C] The resulting expression for the three-body 



loss rate constant can be split in contributions of the 
three-body bound free-space spectrum and continuous 
free-space spectrum: 



^(ry = 0) -> 72^3 — (5 bc 
dr] n\3^o m 



(29) 



The bound-state contribution naturally appears as a 
(rapidly converging) discrete sum over the trimer states: 



bound 



so 



n>0 



(3E t e 



exp 



(pEke 



— 27zn/\sQ 



(30) 

This allows to predict the mean number A^rim of trimers 
with energy -Etrim = — £' t e _27rn /' So ' in the loss-less system 
at thermal equilibrium: Since the contribution to dN / dt 
(to first order in rf) of the term of index n in (130[) is intu- 
itively — 3rtrimAt r i m , where the decay rate of the trimer 
is r trim ^ (2?7/?i|so|)<9i nfl , t £ , t r i m , we obtain 



N nA 3 ^0 



3 3 / 2 (nA 3 ) 2 e-^-" 



(31) 



This agrees with Eq. (188) of [5(| obtained from a chem- 
ical equilibrium reasoning. 

The continuous-spectrum contribution to (|29|) natu- 
rally appears as an integral over positive energies E, see 
Appendix [Cl Mathematically, it can also be turned into 
an easier to evaluate (rapidly converging) discrete sum 



S c 



1 



n>l 



r(i- in \ So \)(pE t y 



(32) 

As expected, S cont is a log-periodic function of E t . In 
practice, due to |so| > 1, it has weak amplitude oscilla- 
tions, between the extreme values ~ 0.478 and ~ 0.522. 
Our continuous-spectrum contribution to L3 is equiva- 
lent, to first order in 77, to the result of a direct three- 
body loss rate calculation for the thermal ensemble of 
free-space three-boson scattering states (56|. 

In experiments, the interaction potential has a finite 
range b, and the actual L3 will deviate from the above 
results. For clarity, we now denote with a star the quan- 
tities corresponding to a finite b. Due to the three-body 
losses, the so-called weakly bound trimer states are actu- 
ally not bound states, they are resonances with complex 
energies E* - ihT*j2. Assuming that T* < \E*\, we 
can name these resonances quasi-bound states or quasi- 
trimers. Their contribution to the decay rate of the Bose 
gas, from the reasoning below Eq. (|3"0)l , can be estimated 
as 



quasi — bound 



3 3 / 2 (nA 3 ) 2 ivE r ^ 



0E* 



(33) 



This is meaningful provided that the thermal equilibrium 
trimer population formula Eq. (188) of (Hoj makes sense 
in presence of losses, that is the formation rate of quasi- 
trimers of quantum number n has to remain much larger 



7 



than r* (in the zero-range framework, this is ensured by 
first taking the limit 77 — > then the limit of vanishing 
density nX 3 —> 0). Evaluation of the finite-6 positive- 
energy continuous spectrum contribution L 



3,cont>0 



to the 

three-body loss rate constant is beyond the scope of this 
work. We can simply point out that, taking the limit 
b — > (with a fixed, infinite scattering length) makes 



3,cont>0 



converge to the value obtained in the zero-range 



finite rj model; further taking the zero-7/ limit gives 



lim — ( lim L 



rj— >0 T] \ b^rO 



7 3,cont>0 



dL 



3,cont 



(V = 0). (34) 



In practice, as soon as b <C X and r\ -C 1, we expect that 



VI. ARBITRARY MIXTURE 

In this Section we consider a mixture of bosonic and/or 
fermionic atoms with an arbitrary number of spin com- 
ponents. The N particles are thus divided into groups, 
each group corresponding to a given chemical species 
and to a given spin state. We label these groups by 
an integer a G {1, . . . , n}. Assuming that there are no 
spin-changing collisions, the number N a of atoms in each 
group is fixed, and one can consider that particle i be- 
longs to the group cr if i € I a , where the I„ 's are a fixed 
partition of {1,...,N} which can be chosen arbitrar- 
ily. For example, a possible choice is I\ = {l,...,Ni}; 
h = {Ni + l,...,iVi + N 2 }; etc. The wavefunction 
ip(rx, ■ ■ ■ , rjy) is then symmetric (resp. antisymmetric) 
with respect to the exchange of two particles belong- 
ing to the same group I a of bosonic (resp. fermionic) 
particles. Each atom has a mass rrii and is subject to 
a trapping potential L^(i"j), and the scattering length 
between atoms i and j is a%j- We set m, = m a and 
a-ij = o-aa' for * £ I<j and j € I a i. The reduced masses 
are = m a m a i j (m a + m a i). We shall denote by P aa i 
the set of all pairs of particles with one particle in group 
a and the other one in group cr', each pair being counted 
only once: 



P aa , = {(*, j) e (la xI a ,)U (/„, xl a )/i< j} . (35) 



The definition of the zero-range model is modified as 
follows: In the contact conditions (|H6p . the scattering 
length a is replaced by cijj, and the limit — > is 
taken for a fixed center of mass position Cy = (mjrj + 
rrijYj) / (rrii + raj); moreover Schrodinger's equation be- 
comes 



N r 



E 



-—A ri +Ui{Ti) 
2m, 



ip = Eip. 



(36) 



troduced the notation in dimension d: 



E 



I] d d r k 



4»*( ClJ ,(r fc ) MlJ )^( Cij, ( r fc)fc#»j)- (37) 

Since a a<T i = a a 'a there are only n(n + l)/2 independent 
scattering lengths, and the partial derivatives with re- 
spect to one of these independent scattering lengths are 
taken while keeping fixed the other independent scat- 



tering lengths. We note that, in Ref. [28j. [Tab. [TTTl 
Eqs. (4a, 4b)] were already partially obtained [7y]. 

In 3D the three-body Efimov effect can occur, e.g. if 
the mixture contains a bosonic group, or at least three 
fermionic groups, or two fermionic groups with a mass ra- 
tio strictly larger than the critical value 13.6069 . . . [57h - 
l59l . In this for single-component bosons, the 

derivatives with respect to any scattering length have 
to be taken for fixed three-body parameter (s), and the 
relation between E and the momentum distribution 
(line 4 of Table HIT]) breaks down, which was not realized 
in [28| When the Efimov effect occurs, there appear 
new relations analogous to the ones given in Section IIVI 
for bosons; we do not give these here for conciseness. 

The derivations of the relations of Tab. Mil are anal- 
ogous to the ones already given for two-component 
fermions and single-component bosons. The lemmas [Ar- 
ticle I, Eqs. (33,35)] are replaced by 



2ttH 2 



,( 2 ) 



(A^\A^) aa , in 3D 



Ha^/a^,)(A^\A^U, in 2D, 



(38) 



where "01 and 1P2 obey the same contact conditions (in- 
cluding the three-body ones if there is an Efimov effect), 
except for the independent scattering length cw , that is 

(i) ■ 

equal to a aa , for ipi, i — 1,2. The momentum distribu- 
tion for the goup a is normalized as J n a (k)d d k/(27r) d = 
N a . The pair distribution function is now defined by 



J2) 



,(u,v) = / d d r 1 ...d d r N \ip(r 1 ,...,v N ) 



x 5(u-Ti)8(y- rj ). (39) 

The Hamiltonian of the lattice model used in some of the 
derivations now reads 



Hiatt = H a+ g ®,<><,> W a 



(40) 



where iJ = Eti -£7^+^^) 



with the discrete 



Our results are summarized in Table [TITl where we in- Laplacian defined by (r|A r |k) = — /c 2 (r|k) (for k in the 
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Three dimensions 



Two dimensions 



dE 



2-kK 1 



d(-l/a 



{A,A) 



(la) 



dE 



irtT 



<9(lna„ 



(lb) 



C a = lim fcV(k) = V(l + <5 a 



,87T/X a 



as 



ft 2 a(-i/<w) 



(2a) 



lim fc 4 n CT (k) = y(l + 5 CTa 



as 



ft 2 <9(lna CTCT /) 



(2b) 



(2) 



R + 

(1 + <5 CT 



-r,R 



dE 1 



v era' 



R + 



27rft 2 9(-l/a CTCT /) r 2 



(3a) 



-r, R 

dE 



m CT + m CT 



to ( - 1 + 5 ^ \ffi a(lna CTCT ,) 



(3b) 



E-E, 



+E 



, = y- 1 dE 

trap ^ q ^ d{-l/a aa ,) 



E 



as 



a(lna CT 



/a^/Ae 7 

m 



f d 3 k 


h 2 k 2 


n CT (k) - 


C CT " 


1 (2tt)3 


2m CT 


" fe4 . 



(4a) 



+ E 



cr<cr' 

d 2 k h 2 k 2 

k<A (2^) 2 2m CT 



n CT (k) 



(4b) 



a 2 £„ 



_ /2vrft : 
2a(-l/<w) 2 ~ V7£ 



2\ 2 



E 



\(A^'\A^) C 

En — E n i 



(5a) 



d 2 E n 



2d{\na aa ,f 



Tih 2 



E 

n',E ,^B„ 



(5b) 



a 2 F 



ac-i/cw) 1 



< 



(6a) 



8 2 F 



a(ln a CT 



< 



(6b) 



d 2 E 



a(-i/a CTCT ,) 5 



< 



(7a) 



d 2 E 



d(\na aa >) 2 y s 



< 



(7b) 



TABLE III: Main results for an arbitrary mixture. In three dimensions, if the Efimov effect occurs, the derivatives must be 
taken for fixed three-body parameter(s), the expression for E in line 4 breaks down, and the last two lines, with derivatives of 
the free energy F and of the mean energy E respectively taken at fixed temperature T and entropy S, are meaningless in the 
absence of spectral selection (see Sec.|V]). 7 = 0.577215 ... is Euler's constant. 



first Brillouin zone) and W aa i — J2(i j)eP ,^ i ,r j b^ d - 
In the formulas of Article I involving the two-body 
scattering problem, one has to replace go by 30,0-0-' , a 
by a aa i and m by 2/Vo'- Denoting the correspond- 
ing zero-energy scattering wavefunction by <j) aa i(T), the 
lemma [Article I, Eq. (56)] is replaced by {%l>'\W aa i\ij)) — 
|<^w(0)| 2 (A',A) aa ,. 

VII. CONCLUSION 

In dimensions two and three, we obtained several re- 
lations valid for any eigenstate of the iV-boson problem 
with zero-range interactions. The interactions are char- 
acterized by the 2D or 3D two-body s-wave scattering 
length a and, in 3D when the Efimov effect takes place, 
by a three-body parameter R t . Our expressions relate 
various observables to derivatives of the energy with re- 
spect to these interaction parameters. Some of the ex- 
pressions, initially obtained in |3lj . were derived in [32| 
with a different technique. For completeness, we have 
also generalized some of the relations to arbitrary mix- 
tures of Bose and/or Fermi gases. 

For the bosons in 3D, especially interesting are the rela- 
tions involving the derivative of the energy with respect 
to the three-body parameter. Physically, one of then 
predicts (to first order in the inelasticity parameter rf) 
the decay rate T of the system due to three-body losses, 



that occur in cold atom experiments by recombination to 
deeply bound dimers. This means that one can extract T 
from the eigenenergies of a purely loss-less (77 = 0) model. 
As an application, we analytically obtained (within the 
zero-range model, and to first order in rf) the three-body 
loss rate constant L3 for the 3D non-degenerate Bose gas 
at infinite scattering length. Experimentally, this quan- 
tity is under current study with real atomic gases [56j. 

Mathematically, the 3D relations hold under the as- 
sumption that the two-body scattering length and the 
three-body parameter are sufficient to make the iV-boson 
problem well-defined, with a self-adjoint Hamiltonian. 
Therefore they may be used to numerically test this as- 
sumption, for example by checking the consistency be- 
tween the values of the derivative of the energy with re- 
spect to the three-body parameter obtained in different 
ways. Three possible ways are: numerical differentiation 
of the energy, the present relation on the short-distance 
triplet distribution function, or the virial theorem which 
also involves this derivative 1601. 
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Appendix A: Derivation of a lemma 

Here we derive the lemma (fTB")) for three bosons in the 
zero-range model. The first step is to exp ress the Hamil- 
tonian in hyperspherical coordinates [l2t |61| : Using the 
value of the Jacobian given below Eq. (p~9|) , 



— 3\/3 / dRR 5 I d 5 n I d 3 C 



5 8 T 1 



h 2 ( f°° f 
= -— 3V3|y dRR 5 J d 5 nA c (R,£l) 

-J d 5 nd 3 CA R (n,C) + dRR 5 J d 3 CA n {R,C) 



(Al) 



where 

Ac(R,n) 



= / d 3 C Uj-Ac^- 



(- 



d 2 



5 d 



dRR 5 { + t;— ) V^2 



\dR 2 RdR 



(A2) 



y>n{^*|^^ ; l (as) 



Tn being a differential operator acting on the hyperangles 
and called Laplacian on the hypersphere. 

The quantity Ar can be computed using the following 
simple lemma: If $i(i?) and ^(R) are functions which 
decay quickly at infinity and have no singularity except 
maybe at R = 0, then 



92 5 9. 



o dRR5 {*\ d R 2 ' RdR 



(A4) 



where Fi(R) = R 2 $>i(R). Expressing the right-hand-side 
of (|A4[) thanks to the boundary condition ([3]) then yields 
the desired result (Tl"3"|) . because the other two contribu- 
tions Ac and .Ao both vanish as we now show. 

The quantity Ac{R, rewritten as | J d 3 CVc ■ 
WVc^ — V^VcV'*) w hh the divergence theorem, is 
zero, since the ip^s are regular functions of C for every 
(i?, JT) except on a set of measure zero. 



It remains to show that 

An{R, C) = for any C and R > 0. (A5) 

We will use the fact that ipi and ip2 satisfy the two- 
body boundary condition (fTJ) with the same a, and apply 
lemma [Article I, Eq. (33)]. More precisely, we will show 
that for any smooth function f(R, C) which vanishes in 
a neighborhood of R = 0, 

J dRR 5 J d 3 C f(R,C) 2 An(R,C) = 0; (A6) 

this clearly implies (|A5)I . To show (|A6|) we note that 

- — 3 V3 / dRR 5 d 3 Cf(R,C) 2 An (R, C) 
2rn J Q J 

h 2 f°° f f 
= -^3%/3 y dRR 5 J d 5 n J d 3 c 

{M)^(m-W^M}, (at) 

which can be rewritten as 

d 3 nd 3 r 2 d 3 r 3 {(M)H(ffo) - [iP*i «■ V> 2 ]} 



2m 



3\/3^ dRR 5 J d 5 n J d 3 C 



(A8) 

The first integral in this expression vanishes, as a conse- 
quence of the lemma [Article I, Eq. (33)]. This lemma is 
indeed applicable to the wavefunctions ftpi: They van- 
ish in a neighborhood of R = (see the discussion in 
Article I), moreover they satisfy the two-body boundary 
condition for the same value of the scattering length a (as 
follows from the fact that R varies quadratically with r 
for small r). The second integral in (| A8|) vanishes as well: 
The contribution of the partial derivatives with respect 
to R vanishes as a consequence of lemma (| A4[) , and the 
contribution of Ac vanishes because the fipi's are regular 
functions of C. 



Appendix B: Relation between T and B for any r/ 

Contrarily to the remaining part of the paper, we as- 
sume here that the inelasticity parameter 77 > 0, so 
that the iV-body wavefunction ip obeys contact condi- 
tions given by Eq.Q and by Eq.Q modified according 
to (f20|) . As a consequence, ip is in general an eigenstate 
of H with a complex energy E — ihT/2, where T is the 
decay rate. If tp is normalized to unity at time then 



r 



d 
"dt 



! )(^0). 



(Bl) 
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This can be transformed using the continuity equation, 
valid for all > 0, and that results as usual from 
Schrodinger's equation: 

9 t |^(X,t)| 2 +div x J = (B2) 



where we collected all the particles coordinates in a single 
vector X = (rx, . . . , rjv) with 3iV components, and where 
we introduced the probability current in M. 3N : 



J = — Im ("0*grad x V')- 



(B3) 



To avoid the singularities that appear in ip for three coin- 
ciding particle positions, we introduce exclusion volumes 
B ijk (e) = {X e R 3N /R ijk < e} for all triplets {i,j,k} 
of particles (of hyperradius iJyjfe) in the integral defining 
1 1 1\) 1 1 2 , taking the limit e — > at the end of the calculation. 
This leads to 



r = -lim / d 3N Xd t (\^(X,t)\ 2 ) 



= -lim V / cP^S-J (B4) 

with the surface element d 3N ~ 1 S oriented towards the 
exterior of -B^fc. Here I e is M. 3N minus the union of all 
Bijk(e); it is thus the set of all the X having all the 
Rijk > e. We use the bosonic symmetry and we single 
out the decay rate due to particles 1, 2 and 3: 



r = - 



N(N - l)(N 
3! 



■ lim 



■3JV-1 



SJ. (B5) 



dB 123 (t) 



The integration domain in Eq. (|B5[) . which is the bound- 
ary of i?i23 (e) , is actually a cylinder in M. 3N , and the coor- 
dinates number 10 to 37V of the surface element d 3N ~ 1 S 
are zero, so that one can keep the contribution to the 
probability current of the first 3 particles only: We can 
thus replace d 3N ~ 1 S ■ J with c? 8 S t • J t , the nine-coordinate 
vectors Jt and d a St coinciding with the first nine coor- 
dinates of J and d 3N ~ 1 S. For fixed r^, . . . , rjy we thus 
have to evaluate 

7 (e) = - d s S t -J t = d 3 r 1 d 3 r 2 d 3 r 3 diy ri , r2 ^3 t , 

JR=e J R>e 

(B6) 

where we used the divergence theorem. We then change 
the integration variables from ri,ra,r3 to C,R, with a 
Jacobian given below Eq. (fl9l) , we use the identity 



y div ri (i/)*grad r . ip — c.c.) = div R (?/>*grad R V> — c.c.) 

i=l 

+ ^divc (V>*grad c V> - c.c.) (B7) 

and we apply backwards the divergence theorem, to ob- 
tain 

7 (e) = -3^3 e 5 J d 3 C J d 5 n ^Im[i>*d R -ip] R = e . (B8) 



Using Eq. Q and Eq.Q modified according to (|20[) . we 
finally obtain 



-N(N- l)(N 



2)^| So |sinh(2ry)||i?|| 2 



(B9) 



with||5|| 2 = J d 3 Cd 3 r 4 ...d 3 r N \B(C,r 4 ,...,r N )\ 2 . In 
the limit ij — > 0, ||i?|| 2 tends to its value in the loss-less 
model and we recover [Tab. HU Eq. (3)] using [Tab. [Til 
Eq. (1)]. 



Appendix C: Free space limit of a virial sum 



Here we derive the free-space limit ([28|) of a sum 
over the internal Efimovian eigenenergies E n (uj) for three 
bosons in a harmonic trap with oscillation frequency w, 
interacting in the zero-range limit with infinite scatte ring 
length. A rewriting of the implicit equation for E n of [44[ 
gives, for n E N: 

lmlnr (i±j^yM ln (^y n ^^ (C1) 

We have introduced the notation E n = E n /(Hu>). Also, 
F(z) is the Gamma function and its logarithm lnr(,z) is 
the usual univalued function with a branch cut on the 
real negative axis. The left-hand side of (|C1|) can be 
shown to be a decreasing function of E n , using relation 
8.362(1) of [H3], so that Eq. (JCTJ) determines E n in a 
unique way. The fact that E t , as given by (f27|) . is the 
free space ground trimer binding energy can be checked 
from (|C1[) by a Stirling expansion for E n — > — oo. 

To evaluate the sum in (|2"5)) for uj — > 0, we collect the 
eigenenergies E n into three groups. The (finite) transi- 
tion group corresponds to |£Vi| not much larger than fno, 
and gives a vanishing contribution to (|28[) for lj — > 0. 
The bound state group corresponds to negative eigenen- 
ergies with |_E„| >• frw; the corresponding free space 
trimer sizes are much smaller than the harmonic oscil- 
lator length [^/(mw)] 1 / 2 , so that the trapping potential 
has a negligible effect and E n (ui) is close to the free space 
trimer energy of quantum number n: 



-E t e 



-27rn/|s | 



(C2) 



This directly leads to the contribution Sbound m (EH])- 

Finally, the third group contains the positive eigenen- 
ergies with E n 3> huj, that shall reconstruct the free space 
continuous spectrum for lj — > 0. As shown in Sec. 3. 3. a 
of [12|, these E n are almost equally spaced by 2Hlo. We 
need here the leading order deviation from equispacing, 
that we parametrize with a "quantum defect" A as 



E n = 2n 

71— > + 00 



A(E n ) + 0{l/n). 



(C3) 



For E n — > +oo, Stirling's formula cannot be immediately 
applied to (jCll) since the argument of the Gamma func- 
tion remains at finite distance from the poles of T (on 
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the real negative axis). Using T(z)T(l — z) = it/ sin(7rz) 
[62| , we obtain the useful identity: 



Im In r ( — ^—^ I = Im In T 



1-sq + E 



+ —E + Im In 



1 



e -T\s \ e -inE 



(C4) 



for all real E. Note that the logarithm in the last term 
of that expression is unambiguously defined (by a series 
expansion of ln(l + it) with u) since e _7r ' s °' < 1. Stirling's 
expansion can now be used in the right-hand side of (|C4[) , 
turning (|Cip into an implicit equation for the "quantum 
defect" A: 



A(E) 



so 



In 



-Im In 



1 



-7r|s L-iirA(E) 



(C5) 

Since exp(— 7r|so|) <C 1, we have a small-deviation prop- 
erty: A(E) only slightly deviates, by 0[exp(— 7r|so|)], 
from the first term in the right-hand side of (|C5I) . This 
deviation was not fully taken into account in §3. 3. a of 



12J. To remain exact, we multiply (IC5[) by in on both 
sides, and we exponentiate the resulting equation. Since 
cxp[— 2ilm ln(l + u)] = (1 + u*)/(l + u), we obtain a solv- 
able equation for exp(i7rA) that determines A modulo 2. 
From the small-deviation property stated above, we can 
lift the modulo 2 uncertainty: 



A(E) = 



In 



-Im In 



-i|sol 



(C6) 

Finally, it remains in (|28[) to replace the sum over n (for 
E n in the third group) by an integral J +o ° dE/(2huj), 
where 2fiw is the leading order level spacing, to obtain 
the continuous spectrum contribution 



cont 



3 3 / 2 
'2k B T. 



dEe 



dA{E) 
d(ln R t 



y d(\nR t ) 

(C7) 

After expansion of b\ n R t A(E) in powers of e 77 ' s ° ' , the 
integral over E can be expressed in terms of the Gamma 
function, which eventually leads to (1321) . 
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